function [ ER ] = ER_partiel( u,sigma,v,tv )

% u et sigma sont les déplacements et contraintes réels
% v et tv sont les déplacements et contraintes virtuels


sigma_v=sigma;
tv_u=tv;

sigma_v(:,2)=sigma(:,2).*v(:,2);
tv_u(:,2)=tv(:,2).*u(:,2);

ER = trapz(sigma_v(:,1) , sigma_v(:,2) - tv_u(:,2) );

integ = 0;

% for i=1:size(u,1)-1
%     
%    integ = integ + ( sigma_v(i+1,2) - tv_u(i+1,2) + sigma_v(i,2) - tv_u(i,2) )*(2*(u(i+1,1)-u(i,1)));
% 
% end
% 
% 
% ER=integ;

% plot(sigma_v(:,1),[sigma_v(:,2),tv_u(:,2)]);
% 
% plot(sigma_v(:,1),sigma_v(:,2)-tv_u(:,2));


end

